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Abstract 

Electrostatics  of  globular  proteins  provides  structural  integrity  as  well  as  specificity  of  biological  function.  This  dual 
role  is  particularly  striking  for  ricin  A-chain  (RTA),  an  A'-glycosidase  which  hydrolyzes  a  single  adenine  base  from  a 
conserved  region  of  rRNA.  The  reported  X-ray  crystallographic  structure  of  the  RTA  mutant  E177A  demonstrated  a 
remarkable  rescue  of  charge  balance  in  the  active  site,  achieved  by  the  rotation  of  a  second  glutamic  acid  (Glu-208) 
into  the  vacated  space.  To  understand  this  conformational  reorganization,  molecular-dynamics  simulations  were 
applied  to  estimate  relative  free  energies  that  govern  the  thermodynamic  stability  of  E177A  together  with  mutants 
E177Q  and  E177D.  The  simulations  anticipate  that  while  E177A  is  a  non-conservative  substitution,  the  protein  is 
more  stable  than  the  other  two  mutants.  However,  the  structural  plasticity  of  the  RTA  active  site  is  not  obtained 
penalty-free,  rather  E177A  among  the  mutants  shows  the  largest  unfavorable  net  change  in  the  electrostatic 
contribution  to  folding.  Of  the  E177A  folded  state,  reorganization  of  Glu-208  lowers  the  electrostatic  cost  of  the 
free-energy  change,  yet  interestingly,  protein  interactions  oppose  the  rotational  shift,  while  solvent  effects  favor  the 
transition.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 

Keywords:  Ricin  A-chain;  Thermodynamic  integration;  Molecular  dynamics;  Hydration 


1.  Introduction 

Knowledge  of  the  determinants  that  govern 
thermodynamic  stability  of  globular  proteins  is 
essential  for  understanding  structure -function 
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correlates.  For  the  large  family  of  ribosome¬ 
inactivating  proteins  (RIPs),  this  includes  factors 
that  control  the  chain-folding  process  underlying 
translocation  across  cell  membranes,  as  well  as 
substrate  recognition  and  binding  to  the  riboso- 
mal  RNA.  Upon  cell  uptake,  RIPs  function  as 
A'-glycosidases  by  removing  an  adenine  base  from 
28S  rRNA,  consequently  leading  to  inhibition  of 
protein  synthesis.  Among  the  RIPs,  ricin  is  the 


0301-4622/01/$  -  see  front  matter  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
PII:  S0301-4622(01)00172-7 


220 


M.A.  Olson  / Biophysical  Chemistry  91  (2001)  219-229 


most  cytotoxic  [1,2].  Ricin  is  a  heterodimer  con¬ 
sisting  of  a  267-residue  catalytic  A-chain  (RTA) 
linked  by  a  disulfide  bond  to  a  262-residue  B- 
chain,  which  mediates  endocytosis  by  binding 
galactose  residues  present  on  various  cell-surface 
glycoproteins  and  glycolipids. 

The  X-ray  ciystal  structures  of  ricin  [3-5]  and 
several  mutants  of  RTA  [6-8]  reveal  that  the 
active-site  cleft  is  balanced  by  a  network  of  nega¬ 
tive  and  positive  charges.  The  core  of  the  RTA 
molecule  is  a  long  connecting  a-helix,  which  ex¬ 
hibits  a  distinct  bend  near  its  C-terminal  end, 
allowing  charged  residues  Glu-177  and  Arg-180  to 
form  an  ion-pair  interaction  at  the  molecular 
surface  of  the  active  site  (see  Fig.  1).  Both  residues 
are  invariant  in  all  known  RIPs  and  are  thought 
to  play  crucial  roles  in  the  catalytic  mechanism 
[9],  For  substrate  recognition  and  binding,  there 
are  seven  additional  highly  conserved  polar 
residues  located  in  the  active  site  and  four  flanking 
solvent-exposed  positively  charged  residues. 

Substitutions  that  alter  the  active-site  charge 
balance  lead  to  structural  rearrangements  of  the 
active  site  plus  decreases  in  solubility  as  a  func¬ 
tion  of  pH  as  compared  to  wild-type  (WT)  RTA 
[8].  One  of  the  more  interesting  structures  arises 
from  the  conversion  of  Glu-177  to  alanine 
(E177A)  [7],  This  mutant  shows  a  second  glutamic 
acid  proximal  to  Glu-177  rotated  into  the  active 
site  to  replace  the  deleted  side  chain,  thus  pre¬ 
sumably  rescuing  the  charge  balance  of  the  WT 
structure.  Catalytic  activity  is  also  remarkably  res¬ 
cued  by  the  E177A  mutant,  showing  an  insignifi¬ 
cant  loss  of  20-fold  [10],  whereas  mutants  E177D 
and  E177Q  exhibited  activity  losses  of  80-  and 
180-fold,  respectively  [10,11],  Replacing  Arg-180 
with  histidine  affects  the  global  stability  of  RTA 
if  the  imidazole  ring  is  deprotonated  [8], 

It  is  well  recognized  that  proteins  in  their  na¬ 
tive  conformational  state  are  in  a  delicate  balance 
of  protein  and  solvent  forces  which  can  be  dis¬ 
rupted  or  stabilized  by  relatively  small  changes  in 
the  side  chains  of  key  amino  acid  residues.  What 
is  not  clearly  understood  is  how  interactions 
among  individual  residues  are  integrated  in  pro¬ 
viding  structural  ‘plasticity’.  To  learn  more  about 
the  factors  that  influence  the  stability  of  the  RTA 


active  site,  molecular-dynamics  (MD)  simulations 
were  carried  out  of  the  folded  and  unfolded  (de¬ 
natured)  structures  of  three  RTA  mutants  de¬ 
rived  from  Glu-177.  Free-energy  simulation  meth¬ 
ods  using  a  finite-difference  thermodynamic  inte¬ 
gration  technique  (reviewed  by  Kollman  [12])  were 
applied  to  estimating  the  relative  folding  stability 
of  mutants  E177Q,  E177D  and  E177A.  Simula¬ 
tions  of  the  folded  states  were  initiated  from  the 
WT  crystal  structure  and  the  unfolded  states  from 
an  extended  chain  conformation.  Free-energy 
differences  were  decomposed  into  contributions 
arising  from  protein  and  solvent  interactions.  Al¬ 
though  free  energies  of  denaturation  have  yet  to 
be  determined  for  RTA  and  the  mutants,  simula¬ 
tions  presented  here  provide  new  insights  into  the 
electrostatics  among  active-site  residues  and 
water  molecules.  This  is  particularly  important  for 
predicting  how  a  protein  thought  a  priori  to  be 
both  non-functional  and  less  stable  can  rescue 
activity,  as  illustrated  by  the  mutation  E177A.  A 
practical  implication  of  establishing  a  predictive 
scheme  for  estimating  the  effects  of  structural 
reorganization  is  the  structure-based  design  of 
protein  mutations  with  modified  rRNA  binding 
specificity.  The  reported  calculations  extend  pre¬ 
vious  theoretical  studies  of  RTA  bound  with  a 
29-mer  RNA  hairpin  and  substrate  analogues 
[13-15],  which  demonstrated  binding  determi¬ 
nants  of  significant  electrostatic  contacts,  several 
of  which  mimic  protein  and  solvent  interactions 
predicted  in  stabilizing  the  native  conformation. 


2.  Theory  and  computations 


Free-energy  differences  for  the  folded  and  un¬ 
folded  state  of  the  enzyme  (denoted  as  E)  were 
calculated  based  on  the  following  thermodynamic 
cycle: 


'fol 


A  Gf  ^  »" 

- ►  Emf 


AGf°l 

^  mu  t 


A  GW1? 

mut 


ACT1 

Efol  ►  E tmf 
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where  A GJ4  and  A G”ut  are  the  free  energies  of 
denaturation  for  the  WT  and  mutant,  respec¬ 
tively,  and  AG“”ft  and  AG“t  are  free  energies  for 
the  unfolded  and  folded  state,  respectively,  com¬ 
puted  via  simulation  in  which  Glu-177  is  con¬ 
verted  to  a  substituting  residue  (Fig.  1).  Substitu¬ 
tions  studied  include  Glu  ->  Gin,  Asp  and  Asp  — > 
X  and  X  ->  Ala,  where  X  denotes  an  intermediate 
residue  designed  to  eliminate  the  negative  charge 
of  Asp  by  replacing  the  carboxylate  group  with  an 
aldehyde  moiety.  This  latter  residue  retains  both 
side-chain  volume  and  hydrogen-bonding  poten¬ 
tial,  and  when  combined  with  the  other  substitu¬ 
tions,  permits  a  systematic  determination  of  the 
structural  conditions  that  allow  charge-balance 
reorganization  to  take  place  in  the  active  site. 


The  relative  free-energy  change  is  obtained  from 
the  following: 

A  AG  =  AG7ut  -  AG,r  =  AG,7  -  AG™  .  (1) 

The  computational  method  employed  in  evalu¬ 
ating  Eq.  (1)  is  based  on  a  finite-difference  ther¬ 
modynamic  integration  algorithm  [12,16]  in  which 
the  respective  free-energy  differences  for  the 
folded  and  unfolded  state  are  calculated  by  using 
the  following  master  equation 

where  the  coupling  parameter,  X,  spans  the  inter¬ 
val  from  0  (representing  the  WT  state)  to  1  (rep- 


Fig.  1.  Schematic  illustration  of  a  ricin  A-chain  Protein,  highlighting  the  ion-pair  interaction  between  residues  Glu-177  and 
Arg-180.  Figure  produced  using  the  MOLSCRIPT  program  [33]. 
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resenting  the  mutant),  and  <  >  denotes  an  ensem¬ 
ble  average  over  a  Boltzmann  sample  of  configu¬ 
rations  governed  by  the  Hamiltonian  H(X).  The 
kinetic  component  is  generally  not  included  in  the 
free  energy  and  the  potential  energy  function 
U(X)  is  sampled.  A  statistical  mechanical  pertur¬ 
bation  theory  [17]  can  be  used  to  compute 
8G(X)/8X  numerically  by  calculating  the  differ¬ 
ence  in  free  energy  AG,  for  a  perturbation  8X 
away  from  X, 


of  the  torsional  term,  <p,  and  y  the  phase  angle; 
Atj  and  Btj  are  constants  depending  on  the  types 
of  atoms  i  and  j;  rtj  is  the  distance  separating 
atoms  i  and  j,  q(  (q  -)  is  the  partial  charge  of  atom 
i  (j),  e  is  the  effective  dielectric  constant;  and  Kt 
is  the  harmonic  restoring  force  constant  acting  on 
atom  i  displaced  A r,  in  boundary  region  R. 
Parameterization  of  Eq.  (4)  as  a  function  of  X  was 
implemented  at  the  level  of  the  force-field 
parameters;  e.g.  the  first  term  is  given  by 


AG,  =  - 


(3) 


where  (3  =  1  /kBT  (Boltzmann  constant  times  ab¬ 
solute  temperature).  Calculations  of  Eq.  (3)  in¬ 
cluded  both  forward  (  +  8X)  and  backward  (-8X) 
perturbations  with  numerical  integration  of  Eq. 
(2)  carried  out  via  a  Gaussian -Legendre  quadra¬ 
ture  [18]. 

MD  simulations  were  carried  out  for  each  X 
point  with  trajectories  generated  for  an  active 
region  composed  of  part  of  the  enzyme  inside  a 
12-A  spherical  boundary  centered  on  the  initial 
position  of  C(3  of  Glu-177  as  derived  from  the 
X-ray  crystal  structure  of  the  WT  (PDB  1RTC). 
The  unfolded  state  was  represented  as  an  ex¬ 
tended  conformation.  Protein  interaction  poten¬ 
tials  were  modeled  using  the  following  form  of 
the  consistent  valance  force  field  [19] 


U( R)=  E  Kb(b-b0f+  E  Ka(e-e0)2 

bonds  angles 

+  E  -^r[l  +  cos(/j(p  -  y)] 

dihedrals 


+  E 

i<i 

1 

+ 


r12 

i) 


E  KjA 


Bu 


r2 

'  l  5 


+  E  — 

sra 
i<j  ‘j 


(4) 


where  R  denotes  the  set  of  all  atomic  coordinates 
specifying  the  system  conformation,  b0  and  0O  are 
the  equilibrium  bond  length  and  angle  parame¬ 
ters,  respectively;  Kb  and  Ka  are  the  bond  and 
bond-angle  stretching  constants,  respectively;  Kd 
is  the  torsional  force  constant;  n  is  the  periodicity 


—  [kXfc(XB)  +  (1  —  X)Xfc(X/<)] 

X{b-  [XbQ(XB) 

+  (1  ~\)b0(\A)]}2  (5) 

where  XA  and  XB  refers  to  the  WT  and  mutant 
state,  respectively. 

Active  regions  were  immersed  in  an  8-10  A 
deformable  layer  of  solvent  water.  Water  was 
represented  by  the  SPC  force-field  potential  [20], 
Ionization  states  for  residues  Asp,  Glu,  Lys,  Arg 
and  His  correspond  to  a  neutral  pH.  The  active 
regions  consisted  of  64  residues  (1063  protein 
atoms)  and  370  water  molecules  for  the  folded 
structure,  and  eight  residues  (126  protein  atoms) 
and  524  water  molecules  for  the  extended  struc¬ 
ture.  Boundary  regions  contained  203  residues 
and  seven  residues  for  the  folded  and  extended 
structures,  respectively.  Amino  acid  residues  lying 
in  unsolvated  outer  shells  were  kept  harmonically 
restrained  at  their  initial  positions  by  employing  a 
force  constant  corresponding  to  a  positional 
mean-square  fluctuation  of  0.52  A2,  determined 
via  an  average  over  all  protein  atoms  in  a  full 
0.5-ns  MD  simulation  of  the  WT  structure.  Bias¬ 
ing  the  search  of  phase  space  was  necessary, 
particularly  for  the  unfolded  form,  where  without 
restraints  the  volume  of  accessible  microstates 
reflected  highly  flexible  peptides,  rather  than  a 
denatured  protein  structure. 

The  simulation  protocol  consisted  of  all  sys¬ 
tems  simulated  at  10  values  of  X  connecting  the 
initial  state  to  the  mutated  state.  The  Xt  simula¬ 
tions  were  initiated  with  400  cycles  of  minimiza¬ 
tion  using  the  steepest  descent  algorithm  plus 
1000  cycles  of  minimization  via  a  conjugate  gradi¬ 
ent  algorithm  followed  by  a  10-ps  equilibration 
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phase.  The  initial  atomic  velocities  were  assigned 
from  a  Gaussian  distribution  corresponding  to  a 
temperature  of  300  K.  Simulations  for  X2,  ^3>— > 
X10  were  carried  out  after  a  10-ps  equilibration 
started  from  the  final  step  of  the  previous  simula¬ 
tion.  Ensemble  averages  were  determined  from 
simulations  of  60  ps  per  X  value  with  non-bonded 
interactions  smoothed  to  zero  beyond  11.0  A  and 
a  constant  dielectric  constant  (e  =  1)  was  used. 
End-point  structures  of  the  mutants  (X  =  1)  were 
determined  following  the  X10  simulation  by  aver¬ 
aging  the  trajectories  of  an  additional  100-ps 
simulation.  Integration  of  the  equations  of  mo¬ 
tion  was  performed  using  a  1.0-fs  timestep. 


3.  Results  and  discussion 

The  free-energy  simulation  results  for  substitu¬ 
tions  Glu-177  — >  Gin,  Asp,  and  Asp  — >  X  and  X 
->  Ala  are  presented  in  Table  1.  Supplementing 
the  calculations  of  AG  are  free-energy  decompo¬ 
sitions  at  the  protein  (AGpro)  and  solvent  (AGsol) 
level,  with  further  decompositions  at  the  force- 


field  level.  Since  we  are  primarily  interested  in 
electrostatic  interactions  of  the  active  site,  van 
der  Waals  and  covalent  interactions  are  com¬ 
bined  into  a  single  non-electrostatic  component. 
Decompositions  derived  directly  from  force-field 
interactions  must  invoke  the  caveat  that  their 
free-energy  contributions  depend  on  the  particu¬ 
lar  integration  path  (see  Boresch  et  al.  [21]  and 
references  cited  therein),  as  well  as  the  choice  of 
force  field.  By  contrast,  the  overall  free-energy 
difference  between  two  end  states  is  independent 
of  the  path  connecting  them.  All  decompositions 
are  consistent  with  Eq.  (2),  with  the  gradient  of 
the  potential  evaluated  via  a  second-order 
finite-difference  scheme. 

The  simulations  show  the  mutants  E177Q, 
E177D  and  E177A  are  destabilized  by  2.0, 1.2  and 
1.1  kcal/mol,  respectively.  Free-energy  differ¬ 
ences  for  E177A  were  calculated  from  adding  the 
contributions  of  E177D,  D177X  and  X177A.  As 
discussed  below,  each  net  change  in  stability  is 
due  to  the  partial  cancellation  of  many  terms, 
both  positive  and  negative,  some  of  which  stabilize 
the  WT  structure  and  others  the  mutants.  These 


Table  1 

Changes  in  free  energies  (kcal/mol)  for  ricin  A-chain  mutants 


Mutation  Free  energy  AG^",f 

contribution  7C  ~7Z  TTTGI 


Total 

Protein 

Solvent 

Total 

Protein 

Solvent 

Total 

Protein 

Solvent 

Total 

Protein 

Solvent 


Total 

Protein 

Solvent 


AG™ 

A  AG 

non-ele 

tot 

ele 

non-ele 

tot 

ele 

non-ele 

6.7 

71.3 

62.2 

9.1 

-2.0 

0.4 

-2.4 

9.5 

47.9 

35.8 

12.1 

-57.3 

-54.6 

-2.7 

-2.7 

23.4 

26.4 

-3.0 

55.3 

55.0 

0.3 

-31.8 

11.8 

43.5 

-31.7 

-1.2 

-1.1 

-0.1 

-36.8 

-21.7 

12.8 

-34.5 

15.3 

17.6 

-2.3 

5.0 

33.5 

30.7 

2.8 

-16.5 

-18.7 

2.2 

-1.9 

101.8 

103.6 

-1.8 

-1.5 

-1.4 

-0.1 

-2.9 

74.7 

76.4 

-1.7 

-61.0 

-59.8 

-1.2 

1.0 

27.1 

27.2 

o 

1 

59.5 

58.4 

1.1 

-12.6 

-15.5 

-2.1 

-13.4 

1.6 

0.8 

0.8 

-15.6 

-9.7 

3.6 

-13.3 

1.2 

3.5 

-2.3 

3.0 

-5.8 

-5.7 

-0.1 

0.4 

-2.7 

3.1 

-46.3 

98.1 

145.0 

-46.9 

-1.1 

-1.7 

0.6 

-55.3 

43.3 

92.8 

-49.5 

-44.5 

-38.7 

-5.8 

9.0 

54.8 

52.2 

2.6 

43.4 

37.0 

6.4 
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findings  are  in  accordance  with  the  work  of  others 
[22-26]  and  would  seem  to  underscore  the  char¬ 
acterization  suggested  earlier  by  Karplus  and  co¬ 
workers  [23],  that  essential  elements  of  the  ther¬ 
modynamics  that  govern  protein  stability  are  ‘hid¬ 
den’  in  the  net  free-energy  change,  either  mea¬ 
sured  or  computed. 

As  with  any  application  of  free-energy  simula¬ 
tion  methods,  there  are  various  approximations 
underlying  the  simulation  model  that  can  intro¬ 
duce  errors  in  the  calculations  [26-28],  These 
include  approximations  due  to  the  force  fields, 
truncation  of  long-range  forces,  sphere  radius  of 
the  active-site  region,  and  computational  sam¬ 
pling  of  phase  space.  Additional  approximation  is 
that  the  unfolded  state  lacks  any  residual  sec¬ 
ondary  structure,  although  several  experimental 
and  theoretical  works  may  disfavor  this  assump¬ 
tion.  Requisite  testing  of  the  cutoff  distance  and 
sphere  size  used  in  the  models,  as  well  as  the 
length  of  700  ps  for  the  MD  simulations,  suggest 
that  they  were  adequate.  Systematic  errors  were 
estimated  from  analyzing  the  hysteresis  derived 
via  forward  and  backward  perturbations,  and 
showed  values  of  approximately  +0.6  kcal/mol 
for  each  net  A G  determination. 

3.1.  Mutation  Glu-177  ->  Gin 

The  individual  free-energy  changes  from  the 
amide  substitution  reflect  predominately  elec¬ 
trostatic  components  and  show  large  unfavorable 
contributions,  with  values  of  62.6  kcal/mol  for 
the  unfolded  state  and  62.2  kcal/mol  for  the 
folded  state.  The  net  effect  on  A  AG  from  elec¬ 
trostatic  interactions  is  to  marginally  stabilize 
E177Q  by  0.4  kcal/mol.  Destabilization  arises 
from  non-electrostatic  terms  for  a  loss  of  -2.4 
kcal/mol.  Of  the  total  A  AG,  electrostatic  inter¬ 
actions  from  the  protein  are  found  to  stabilize  the 
unfolded  chain  conformation  and  destabilize  the 
folded  form,  yielding  a  net  decrease  in  stability  of 
-54.6  kcal/mol.  This  loss  is  offset  by  a  favorable 
electrostatic  desolvation  free  energy  of  glutamine 
vs.  glutamate,  contributing  55.0  kcal/mol  toward 
AAG. 

The  average  MD  simulation  results  for  the 
active  site  of  the  WT  structure  and  E177Q  are 


illustrated  in  Fig.  2a, b.  Overall,  the  modeled  WT 
structure  is  in  good  agreement  with  the  crystallo¬ 
graphic  determination  reported  by  Robertus  and 
coworkers  [3],  The  carboxylate  of  Glu-177  is  an¬ 
chored  in  place  by  an  ion-pair  interaction  between 
Oe2  and  Arg-180  NH2,  and  a  strong  hydrogen 
bond  between  Osl  and  the  hydroxyl  of  Tyr-123. 
The  charge  neutralization  arising  from  E177Q 
alters  both  interactions.  Moreover,  the  perturba¬ 
tion  leads  to  unfavorable  changes  in  hydration 
surrounding  the  active  site.  Three  hydrogen  bonds 
with  water  are  calculated  for  the  WT  Glu-177; 
Oel  bonds  to  two  bridge  forming  waters,  one 
forming  a  bridge  to  NH2  of  Arg-180  (similar  to 
the  water  denoted  as  no.  334  of  the  crystal  struc¬ 
ture),  and  a  second  with  Oe2  of  Glu-208;  Oe2  is 
hydrogen  bonded  to  an  additional  third  water. 
The  simulation  structure  of  E177Q  shows  one 
weak  binding  water  molecule  interacting  with  the 
side  chain  of  Gln-177. 

3.2.  Mutation  Glu-177  -* Asp 

Eliminating  one  methylene  group  from  the  side 
chain  of  Glu-177  yields  a  mutant  E177D  that  is 
less  stable  by  1.2  kcal/mol  than  the  WT  struc¬ 
ture,  but  0.8  kcal/mol  more  stable  than  E177Q. 
In  contrast  with  E177Q,  the  net  electrostatic  ef¬ 
fect  of  the  mutation  is  to  favor  the  unfolded  state 
(AAG  =—1.1  kcal/mol).  A  further  difference 
with  E177Q  is  that  the  net  protein  electrostatic 
contribution  for  E177D  favors  the  mutant  form, 
while  solvent  effects  destabilize  the  substitution. 
Of  the  folded  conformation,  the  electrostatic  per¬ 
turbation  on  protein  interactions  is  much  less 
than  the  amide  conversion  (~  13  vs.  36  kcal/mol). 
Fig.  2c  shows  the  carboxylate  of  Asp-177  preserv¬ 
ing  the  ion-pair  interaction  with  NH2  of  Arg-180, 
although  at  an  energetic  cost  of  shifting  the  posi¬ 
tively  charged  residue.  Similarly,  a  hydrogen  bond 
is  formed  comparable  to  the  WT  structure  with 
the  hydroxyl  of  Tyr-123,  yet/equires  a  conforma¬ 
tional  displacement  of  ~  1  A.  Because  the  substi¬ 
tution  conserves  the  negative  charge,  the  carboxy¬ 
late  of  Glu-208  is  predicted  to  occupy  a  position 
similar  to  that  of  the  WT  structure,  maintaining 
the  interaction  with  Gln-173. 

The  simulation  structure  of  the  active-site  re- 
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Fig.  2.  Results  of  the  molecular  dynamics  simulations  of  wild-type  structure  and  amino  acid  substitutions  at  position  177  of  the  ricin 
A-chain  protein  in  explicit  solvent  water,  (a)  The  active  site  of  the  wild-type  structure.  Water  molecules  are  illustrated  as 
ball-and-stick  models.  Only  the  water  interacting  with  side  chains  of  Glu-177  and  Glu-208  are  shown,  (b)  Contrast  between  the 
calculated  E177Q  active  site  and  the  average  simulation  wild-type  structure.  The  thick  line  denotes  the  mutant  and  thin  line  shows 
the  wild-type  protein.  Waters  correspond  to  the  mutant  structure.  Simulation  structures  for  (c)  E177D  and  (d)  E177A  mutants. 
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gion  indicates  two  waters  are  bound  by  Asp- 177 
and,  because  of  the  structural  reorganization  in 
residues  Tyr-123  and  Arg-180,  the  solvent  con¬ 
tribution  is  less  favorable  than  that  calculated  for 
E177Q.  Nevertheless,  the  bound  waters  are  more 
consistent  with  the  WT  structure  and  should  con¬ 
tribute  favorably  to  catalysis.  For  the  unfolded 
state,  displacement  of  the  charge  yields  a  more 
stable  free  energy  of  solvation  than  that  corre¬ 
sponding  to  the  charge  deletion  of  E177Q. 

3.3.  Mutation  Glu-177  -» Ala 

The  free-energy  perturbation  of  glutamate  to 
alanine  was  expanded  into  three  intermediate 
steps:  E177D;  D177X;  and  X177A.  As  mentioned 
above,  X  denotes  a  side  chain  in  which  the  car- 
boxylate  of  Asp  has  been  replaced  with  an  alde¬ 
hyde  group.  The  idea  was  to  design  a  residue  that 
eliminates  the  forking  interactions  of  the  carboxy- 
late  while  retaining  the  potential  for  a  single 
hydrogen  bond.  From  Table  1,  the  conversion  of 
Asp  to  X  significantly  disrupts  the  electrostatic 
stability  (A A G  =  - 1.4  kcal/mol)  greater  than  ei¬ 
ther  E177Q  (  +  0.4  kcal/mol)  or  E177D  (-1.1 
kcal/mol).  Protein  contributions  are  unfavorable 
in  both  unfolded  and  folded  conformations,  lead¬ 
ing  to  a  net  destabilization  of  -61.0  kcal/mol. 
The  calculations  for  the  folded  form  show  that 
the  decreased  polarity  for  hydrogen  bonding  of  X 


contributes  to  an  electrostatic  free-energy  loss  of 
nearly  30  kcal/mol  greater  than  the  charge  neu¬ 
tralization  of  E177Q,  and  approximately  60 
kcal/mol  for  a  comparison  with  E177D. 

Structurally,  the  residue  X-177  preserves  the 
hydrogen  bond  with  Arg-180  (data  not  shown), 
although  much  weaker  than  the  WT  ion  pair.  The 
substitution  also  disrupts  the  electrostatic  attrac¬ 
tion  between  X-177  and  Tyr-123.  For  residue 
Glu-208,  the  calculated  structure  indicates  a  lack 
of  significant  dihedral  angle  perturbation.  This 
implies  that  the  reduced  solvent-exposed  polarity 
at  position  177  while  retaining  side-chain  volume 
is  insufficient  to  allow  the  proximal  glutamate  to 
reorganize  the  charge  balance  of  the  active  site. 
As  for  solvent  effects,  the  net  contribution  of 
transferring  the  mutant  from  vacuum  to  solution 
is  favorable  ( —  60  kcal / mol)  and  is  similar  to 
E177Q  (~  55  kcal/mol). 

The  replacement  of  X  with  alanine  is  the  final 
step  in  obtaining  E177A  and  yields  a  protein 
X177A  that  shows  an  increase  in  thermodynamic 
stability  by  a  A  AG  of  1.6  kcal/mol.  The  signifi¬ 
cance  of  this  result  is  that  the  improved  stability 
is  obtained  by  the  dihedral  rotation  of  Glu-208  to 
fill  the  active-site  cavity.  The  recovery  in  stability 
yields  the  least  unfavorable  electrostatic  compo¬ 
nent  among  the  series  of  folded  protein  changes 
arising  from  the  substitutions.  Moreover,  the 
modification  displays  favorable  net  protein  and 
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Fig.  3.  Molecular  modeling  drawing  of  the  reorganization 
calculated  for  the  E177A  structure  and  its  comparison  with 
the  wild-type  structure  for  ricin  A-chain.  The  color  yellow 
designates  the  side  chains  of  E177A  and  the  blue  colored 
worm  shows  the  backbone.  Side  chains  of  the  wild-type  struc¬ 
ture  are  shown  in  white,  with  oxygen  atoms  illustrated  in  red 
and  nitrogen  atoms  in  blue.  Figure  produced  using  the  GRASP 
program  [34]. 

solvent  contributions  for  both  unfolded  and  folded 
conformations,  however,  the  total  electrostatic 
component  from  solvent  interactions  is  unstable. 

For  E177A,  individual  contributions  to  A  AG 
are  dominated  by  unfavorable  electrostatic 
changes,  yet  interestingly,  it  is  the  non-electro- 
static  component  that  makes  E177A  more  stable 
than  either  E177Q  or  E177D.  The  average  simu¬ 
lation  structure  for  the  active  site  of  mutant 
E177A  is  shown  in  Fig.  2d  with  a  superposition  of 
WT  simulation  structure  and  the  two  structures 
are  contrasted  schematically  in  Fig.  3.  The  struc¬ 
tural  difference  is  primarily  the  rotation  of  the 
carboxylate  of  Glu-208  in  accordance  with  the 
X-ray  crystal  structure  [7],  As  expected,  Ala-177 
has  no  strongly  interacting  waters  in  the  folded 
state  and  only  one  in  the  unfolded  state  through 
backbone  interactions.  Presumably,  the  functional 
hydration  sites  for  catalysis  are  now  provided  by 
the  side-chain  reorganization  in  Glu-208. 

From  the  average  E177A  simulation  structure, 


the  side  chain  of  Glu-208  has  dihedral  angles  of 
X,  =  - 140.2°  and  =  62.5°,  which  correspond  to 
a  less-populated  rotameric  state  [29],  Because  of 
the  rotation,  the  carboxylate  of  208  breaks  the 
WT  interaction  with  Gln-173,  yet  gains  two  alter¬ 
native  hydrogen  bonds  within  the  protein  active 
site.  The  first  is  with  Tyr-123,  which  undergoes  a 
rearrangement  allowing  the  interaction  to  take 
place,  and  the  second  is  with  the  guanidinum 
group  of  Arg-180,  which  appears  to  move  upward 
in  the  active  site.  The  X-ray  crystal  structure  of 
the  mutant  E177A  also  shows  Tyr-80  rotated  out¬ 
ward  toward  the  solvent  [7],  whereas  the  simula¬ 
tion  of  700  ps  places  the  residue  in  a  similar 
conformation  as  found  in  the  WT  structure.  Owing 
to  this  lack  of  significant  movement  in  Tyr-80  in 
the  simulation,  it  is  anticipated  that  the  calcu¬ 
lated  A  AG  of  -1.1  kcal/mol  may  approximate 
an  upper  bound  in  stability  change  for  E177A. 

To  understand  better  the  role  of  Glu-208  and 
its  rotation  into  the  active  site,  additional  free- 
energy  simulations  of  X177  ->  A  were  preformed 
where  the  glutamate  side  chain  was  restricted  to 
its  WT  configuration.  The  motivation  for  these 
calculations  was  a  crystallographic  analysis  of  the 
E177A  structure,  which  showed  that  approxi¬ 
mately  15%  of  the  molecules  retained  the  non- 
rotated  configuration  for  Glu-208  [7].  Using  occu¬ 
pancy  populations  of  85  and  15%,  a  free-energy 
difference  between  the  two  Glu-208  configura¬ 
tions  is  estimated  at - 1  kcal/mol,  favoring 

side-chain  rotation.  From  the  observation  that  the 
side-chain  rotation  of  Glu-208  occurs  during  the 
X  — >  A  substitution,  the  simulation  model  for  this 
perturbation  provides  a  reasonable  approach  for 
calculating  the  conformational  free-energy  dif¬ 
ference.  Simulation  results  for  the  reorganization 
are  presented  in  Table  2.  Because  of  energy  con¬ 
servation  problems  with  the  WT  configuration, 
thermal  initialization  was  extended  for  both  states 
to  75-80  ps  per  X  point  starting  with  the  \7 
perturbation. 

The  results  show  that  the  side-chain  reorgani¬ 
zation  in  Glu-208  is  favored  thermodynamically 
by  a  AG  of  -1.4  kcal/mol.  Without  forcing  the 
rotation  in  Tyr-80,  the  calculated  AG  is  in  excel¬ 
lent  agreement  with  the  crystallographic  estimate. 
The  calculations  further  show  that  the  elec- 
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Table  2 

Conformational  reorganization  free  energy  (kcal/mol)  for 
residue  Glu-208  in  the  structure  E177A 


Free  energy 
contribution 

AG 

expa 

tot 

ele 

non-ele 

Total 

-1.0 

-1.4 

-0.5 

-0.9 

Protein 

- 

6.9 

7.0 

-0.1 

Solvent 

- 

-8.3 

-7.5 

-0.8 

“Experimental  estimate  taken  from  Kim  et  al.  (1992). 


trostatic  component  shifts  the  equilibrium 
between  the  two  conformers  toward  structural 
reorganization.  Moreover,  and  perhaps  unex¬ 
pected,  the  results  indicate  that  relative  to  the 
WT  configuration  the  protein  interactions  of  the 
active  site  oppose  the  rotation  of  Glu-208.  De¬ 
spite  the  favorable  pair-wise  interactions  of  Glu- 
208  with  side  chains  of  Arg-180  and  Tyr-123, 
electrostatic  clashes  with  the  residual  protein 
backbone  make  AGpr0  less  than  optimal.  This  is 
not  to  say  that  the  reorganized  protein  free-en- 
ergy  surface  is  unfavorable,  but  rather  the  more 
populated  E177A  active-site  configuration  con¬ 
tains  electrostatic  strain.  Conversely,  solvent  in¬ 
teraction  promotes  the  conformational  reorgani¬ 
zation.  The  folding  of  E177A  retains  greater  sol¬ 
vent  accessibility  for  the  side  chain  of  Glu-208 
than  the  WT  configuration,  and  thus  yields  a 
lower  desolvation  cost  reflected  in  the  calculated 
electrostatic  component  of  AGsol.  A  similar  argu¬ 
ment  can  be  suggested  for  Tyr-80  in  the  crystal 
structure,  where  side-chain  rotation  is  favored 
due  to  solvent  interactions. 

From  the  model  calculations,  a  general  picture 
of  the  thermodynamic  stability  of  the  RTA  active 
site  emerges.  While  the  electrostatic  cost  of  des¬ 
olvation  due  to  the  release  of  water  from  the 
protein  active-site  surface  upon  folding  favors 
charge  deletion  of  Glu-177,  maintaining  ion-pair 
complementarity  is  a  necessary  component  for 
minimizing  the  thermodynamic  loss  in  stability. 
This  is  evident  from  eliminating  the  negative 
charge  at  position  177  without  allowing  significant 
electrostatic  reorganization  (illustrated  by  mu¬ 
tants  E177X  or  the  WT  configuration  of  E177A, 


both  showing  a  value  for  A  AG  <  -2.5  kcal/mol). 
Overall,  structural  plasticity  arises  from  searching 
the  free-energy  surface  for  the  configuration  that 
yields  the  minimum  reorganization  penalty.  Pro¬ 
tein  reorganization  of  RTA  from  charge  deletions 
is  unfavorable  in  terms  of  electrostatics,  while 
solvent  reorganization  shifts  the  populations  on 
the  energy  landscape  in  a  fashion  that  recovers  a 
significant  fraction  of  the  loss  in  conformational 
stability.  Although  still  debated  is  the  issue  of 
whether  hydrogen  bonds  make  stabilizing  con¬ 
tributions  to  the  free-energy  difference  between 
the  unfolded  and  folded  state  [30-32],  it  is  clear 
that  charge  balance  is  a  critical  element  of  op¬ 
timizing  the  self-organization  of  the  native  ter¬ 
tiary  fold. 

The  calculations  presented  here  demonstrate 
that  free-energy  simulation  methods  are  capable 
of  detecting  a  remarkable  reorganization  of  the 
RTA  active  site  subjected  to  the  non-conservative 
perturbation  of  E177  -» A.  The  E177A  crystal 
structure  emphasizes  the  importance  of  consider¬ 
ing  structural  plasticity  in  interpreting  mutagenic 
results.  A  promising  aspect  of  the  calculations  is  a 
modeling  approach  for  understanding  the  under¬ 
lying  forces  that  stabilize  RTA  and  its  mutants. 
Moreover,  the  calculations  provide  a  strategy  for 
anticipating  the  effects  of  additional  substitutions 
aimed  at  eliminating  the  cytotoxicity  of  ricin,  while 
preserving  antigenic  determinants  characteristic 
of  the  WT  protein. 
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